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^^ Abstract 

^^ The mobility of electrons and holes in silicon depends on many parameters. Two of them are the electric field and the 
^^ temperature. It has been observed previously that the mobility in the transition region between ohmic transport and 

saturation velocities is a function of the orientation of the crystal lattice. 

This paper presents a new set of parameters for the mobility as function of temperature and electric field for <111> 

and <100> crystal orientation. These parameters are derived from time of fiight measurements of drifting charge 

carriers in planar p+nn+ diodes in the temperature range between -30°C and 50°C and electric fields of 2x10^ V/cm to 

2x10"^ V/cm. 
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■1. Introduction 
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The transport properties in bulk silicon have been stud- 
ied in great detail [l| . Precise knowledge of the transport 
properties is important, as in applications with large drift 
|>^distances, e.g. radiation sensors, even small uncertainties 
^I^ 'in the drift velocities can have noticeable effects. 
^ Mh In a radiation sensor ionizing radiation creates electron 
hole pairs in the silicon bulk material, which drift in the 
CS| electric field created by an applied voltage and induce a 



> 

(N 

m 

o 
o 



X 
S3 



signal in the readout electrodeqj. Diffusion processes occur 
simultaneously causing a lateral spread of the charge car- 
riers and a drift velocity distribution where some charge 
carriers are accelerated and some are decelerated due to 
|the diffusion. The current induced by each charge carrier 
■is proportional to its drift velocity, which, if diffusive pro- 
cesses are neglected, is given by the product of mobility 
and electric field. Thus, in order to optimize radiation 
sensors, precise knowledge of the mobility as function of 
the electric field is necessary. 

Common radiation sensors are built in such a way that 
drift fields are in the order of 10^ V/cm, making the charge 
collection times of these sensors sensitive to the charge car- 
rier mobility in the transition regime between ohmic and 
saturation behavior. Especially for sensors damaged by 
irradiation the charge collection time is an important op- 
timization parameters as charge carriers are trapped and 
thus degrade the obtained signal. 
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Parameterizations of charge carrier mobilities can be 
found in the literature, however the most recent measure- 
ments using the time of fiight technique have been per- 
formed by Canali et al. in 1971 [2]. There have been 
quite a few reviews of this data, the one by Jacoboni et 
al. [3] is especially useful, as it provides a temperature 
depended parameterization of the drift velocity in <111> 
direction. There have been other reviews incorporating 
other measurements and proposing different parameter- 
izations [4, 5, IgI, I?!, however either temperature depen- 
dence or anisotropy is neglected. For electric field values 
of 10^ V/cm in <111> orientated silicon at a tempera- 
ture of 300 K the aforementioned parameterizations pro- 
duce values scattered by approximately 10% (best case), 
however the scattering increases significantly for different 
temperatures (57% for holes at 250 K). None of the men- 
tioned reviews provides a parameterization for the <100> 
direction. 

In this work measurements of pulse shapes of p+nn+ 
diodes in the temperature range between -30° C and 50° C 
with different applied bias voltages are compared to simu- 
lations in order to determine a set of parameters, describ- 
ing the mobility as function of electric field and tempera- 
ture for both <111> and <100> orientated silicon. 



2. Time of flight method 

The time of fiight method (sometimes also called tran- 
sient current technique, TCT) is a well established tech- 
nique to determine effective trapping times in irradiated 
silicon diodes [sl, |9|. It can also be used to investigate 
charge multiplication or effects of high charge carrier den- 
sities in segmented devices [lO|, lll|. Electron hole pairs 
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are created inside a reverse biased diode either by ionizing 
radiation or laser light. These charge carriers then drift un- 
der the influence of the electric field and induce a current 
in the readout circuit. The current pulse contains infor- 
mation about the amount of drifting charge and their drift 
velocity. By using a laser source with appropriate wave- 
length the expected current pulses for x-rays, heavy ions or 
minimum ionizing particles can be simulated, allowing the 
optimization of application specific readout electronics. 

Planar p+nn+ silicon diodes with <111> and <100> 
crystal orientation were investigated using the time of fiight 
method. 

In order to measure the mobility of electrons and holes 
separately we create electron hole pairs by injecting laser 
light either from the p+ or from the n+ side. The pene- 
tration depth of the light is short compared to the thick- 
ness of the investigated samples, thus the other type of 
charge carriers is collected very quickly (within approxi- 
mately 100 ps). 

To evaluate whether effects of high charge carrier den- 
sities (also called plasma effects) were present the diameter 
of the laser spot on the sample was changed by focusing 
and defocussing. As no changes in the recorded current 
pulses were observed, high density effects can be consid- 
ered negligible. 

3. Experimental setup 

Charge carriers were created with a laser of 660 nm 
wavelength and the time resolved current pulses of the in- 
vestigated diode was read out by a Miteq AM- 1309 wide- 
band amplifier and a Tektronix DPO 7254 2.5 GHz oscil- 
loscope. 

The sample was mounted on a substrate that allowed 
light injection from both sides while providing a stable (± 
0.1 K, rms) temperature in the range of 240 K to 340 K 
and applying the high voltage to the rear side of the diode. 
The guard ring was grounded. 

3.1. Equivalent circuit for SPICE simulations 

An equivalent circuit of the setup was used for SPICE 
simulations (see Figure [2j). In the equivalent circuit de- 
tector capacitance (9.2 pF or 9.4 pF), signal cable length 
(3.05 m), bandwidth of the amplifier (1 GHz) and imper- 
fections of the setup (inductances and capacitances) were 
taken into account. The transfer function of the whole 
setup has been checked and was found to introduce distor- 
tions, thus all simulations have been convoluted with the 
transfer function (shown in Figure [3]). 

3.2. Laser properties 

The laser system was manufactured by Picoquant GmbH 

12| and provided short and intense light pulses with a 

FWHM < 100 ps. The time resolved pulse structure was 

provided by the manufacturer. For this study laser light 

with a wavelength of 660 nm (± 2 nm) was used. The 
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Figure 1: Sketch of the experimental setup showing the polarity of 
the applied voltage. The laser light can be shined in either from the 
junction side or opposite to the junction (indicated by the arrows 
with broken lines). The diode has a thickness of 280 |Ltm or 272 |j.m 
and is connected to the readout circuit. 



pulse energy of 140 pJ was attenuated with an optical at- 
tenuator, which has no effect on the time structure of the 
pulses. 660 nm light has an attenuation length of roughly 
3 |Ltm in silicon at 20° C. 

3.3. Investigated samples 

The investigated samples were planar p+nn+ diodes 
fabricated by CiS Forschungsinstitut fiir Mikrosensorik und 
Photovoltaik GmbH 



13|. The silicon used is high resis- 



tivity n-type diffusion oxygenated float zone silicon with 
<100> and <111> orientation manufactured by Wacker. 
The effective doping of the <100> and <111> samples are 
8.2x10^^ cm~^ and 1.1x10^^ cm~^. The resulting depletion 
voltages are 49 V and 63 V with dark currents at the deple- 
tion voltage of 0.72 nA and 1.64 nA, respectively. The p+ 
implantation has a peak value of 3x10^^ cm~^ and a junc- 
tion depth of 0.5 iitm. The n+ implantation has a peak 
value of 2x10^^ cm~^ and depth of 1 jam. The samples 
have very low concentrations of lattice defects, leading to 
trapping times in the order of milliseconds, which is many 
orders of magnitude larger than the longest pulse dura- 
tion recorded in this work (approximately 35 ns). Thus 
any trapping effects have been considered negligible. 

From the measured pad capacitance of 9.23 pF and the 
pad area of (4.95 mm)^ we obtain a sample thickness of 
280 |LLm using the standard formula for a parallel plate ca- 
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Figure 3: Transfer function derived from the SPICE simulations for 
9.2 pF capacitance. A similar transfer function was calculated for 
9.4 pF capacitance, however the deviations are to small to be seen 
on this scale. The inset shows a zoom of the peak structure of the 
transfer function using the same units on the x- and y-axis. 



Figure 2: Equivalent circuit for SPICE simulations. Dominant ele- 
ments are the diode capacitance (only 9.2 pF shown) in parallel to 
the current source (marked by the rectangle), the system inductance 
(30 nH), the lossy transmission line (3.05 m RG-174 cable) and the 
band pass filtering properties of the amplifier. 



pacitor without edge effects. The corresponding numbers 
for the <111> diode are 9.42 pF, (4.93 mm)^ and 272 M,m. 
For the measurement of the pad capacitance the capaci- 
tance of the guard ring to the backplane and the sensor 
edge is subtracted by the zero adjustment of the capac- 
itance bridge. The remaining effect of the guard ring is 
estimated to be well below 1 %. The estimated uncer- 
tainty of the diode thickness is ± 2 jam. 

Mechanical measurements of the thickness yield some- 
what higher values but include 'dead' layers like passiva- 
tion and aluminizations. The accuracy of the final capac- 
itance has the biggest influence on the determination of 
the thickness, as a change of 0.1 pF changes the measured 
thickness by approximately 3 iitm. 

The samples feature an opening on the top side and a 
grid like aluminum layer on the rear side to allow light in- 
jection. The gap between the metalization of diode guard 
ring is 20 |Ltm wide, the distance between the corresponding 
implantations is 10 iitm. A sketch of the sensor is shown 
in Figure [H 

4. Measurements 

Measurements have been performed for front side illu- 
mination (electron signal) and rear side illumination (hole 
signal). Bias voltages of 100 V, 200 V and 500 V were 
applied and these measurements were done in tempera- 
ture steps of 10 K in the range between -30 °C and 50 °C 



(243.15 K and 323.15 K). A set of measurements and cor- 
responding simulations are shown in Figure IH 

To reduce the random noise component of the system 
averages 4096 pulses were taken. The remaining noise is 
no longer random and has an rms value of < 0.4 |llA. 

5. Simulations 

The simulations were done in two steps. At first the ex- 
pected current pulse was calculated using a custom written 
transport code, then the resulting current pulse was con- 
voluted with the time structure of the light pulse and with 
the transfer function of the setup. 

The custom transport code is taking charge drift and 
diffusion into account. It is a ID calculation on a grid with 
100 nm cell size {Ax) with 5 ps time step (At). 

Above full depletion the electric field inside the diode 
is linear with its highest value at the p+n-j unction {x = 0) 
and described by the following equation: 



E{x) = - ( Ubias 
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dep 



2x 



(1) 



Where x is the position inside the sensor (center of 
corresponding cell), d the sensor thickness (280 iitm or 
272 |Ltm), Uhias the applied bias voltage and Udep the de- 
pletion voltage of the diode. 

The normalized charge carrier distribution {Ne^h{x^ t = 
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Ax 



exp{—x*/Xabs)) generated by the laser light fol- 



lows an exponential function with attenuation length Xabs = 
3 |LLm with x^ = X for front side injection and x* = d — x 
for rear side injection. 

The induced current is calculated according to the Ramo 
theorem [14]. The weighting field of a planar diode is lin- 
ear and the induced current (normalized to the elementary 
charge eo) can be approximated as: 





/xO [cmVVs] 


u^"* [cm/s] 


o<100> 


^<iii> 


electrons 
holes 


1440±15 
474±10 


(1.054±0.038)xl0^ 
(0.940±0.027)xl0^ 


0.992±0.004 
1.181±0.003 


1.075±0.004 
0.924±0.002 



Table 1: Parameters for the field dependent mobility at 300 K. 



Temperature parameter ai for 



electrons 
holes 



-2.260±0.007 
-2.619±0.007 



-0.602±0.003 
-0.226±0.002 



0.572±0.003 
0.644±0.003 



0.220±0.001 
0.550±0.002 



Table 2: Temperature parameters a^ for the field dependent mobility reduction. 



/(t + At) = ^[/e(x,t + At)-4(x,t + At)] (2) 
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5.1. Drift and diffusion 

The charge drift and diffusion was simulated using the 
following basic equations: 



5x{x) = ±iie,h{x)E{x)At 
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Where Sx{x) is the distance traveled within one time 
step (At). It is negative for holes as they drift into the op- 
posite direction of electrons. Average values (e.g. jr^{x)) 
have been calculated by adding the values for the trans- 
versed cell and dividing by the number of cells traveled 

(±l/(5x* YTx ^ •> ^^^ ^^* ^^ ^^ rounded to nearest cell cen- 
ter). 5x was determined by adding the drift times for each 
cell until the total drift time {t drift) exceeded the time 
step when Uceii cells have been transversed, then calcu- 
lating 5x = ^ceii^Ax ^^^ Diffusion is simulated using the 

^ tdrift ^ 

Einstein relation (I^e,/i = ~P^{x) — ). 

The number of charge carriers at each grid cell is cal- 
culated as the integral of a Gaussian function shifted by 
5x with G = ^/2DAt from x - ^ to x + ^ (beginning of 
cell to end of cell). 

These shifted Gaussian functions are calculated for each 
grid cell and summed over all grid cells. 



In order to avoid problems with charge carriers leaving 
the diode (x < or x > d) and diffusing back into the 
diode, a charge carrier which reaches the last grid point 
along the drift direction is removed from the simulation. 

The numerical stability of the simulation has been checked 
by varying the grid size and time step by a factor of up to 
five. No effect on the pulse shape was observed. Charge 
conservation was checked and no charge loss was found. 

5.2. Mobility model 

The field dependence of the mobility was parameter- 
ized in the following way: 



l^e,h{x) 



11^ 
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Where jiie,h{x) is the mobility of electrons or holes for 



a given grid cell x with the electric field E{x)^ /i^ ^ 
ohmic mobility (also called low field mobility), 



^e,h 



the 
the 



saturation velocity and Pe,h a parameter describing the 
transition between ohmic behavior and saturation velocity. 
The temperature dependence of these parameters was 
modeled using a simple exponential law: 



Parameteri{T) = P ar ameteri {300 K) 

V 300iv / 



T V 



(12) 



6. Results 



To determine the free parameters of the model a mod- 
ified x^ procedure is used. At first current pulses for all 
measurements were simulated (<111> and <100> sample 
for 100 V, 200 V and 500 V at ah temperature steps; 54 
pulses in total) and the squares of the difference between 
measurement and simulation were summed up. This x^ 
sum was calculated for every set of parameters and the 
parameters with the smallest x^ are presented in Tables [1] 
and m Errors have been calculated by using the Minuit 
fitting program [15.] of the ROOT package [16]. Transients 
calculated with this set of parameters are shown in Figure 
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Figure 4: Comparison of simulated pulses and measurements at T=30°C. The upper row shows pulses for <111> orientation, the lower for 
<100> orientation; The left column electrons, the right column holes. The measured pulses (solid lines) are generally well reproduced by the 
simulations (dotted: using the mobility parameters of Jacoboni et al. Q, dashed: using the mobility parameters of this work), although some 
small ringing features are not correctly reproduced by the SPICE simulation. For electrons the simulated pulses of both models differ only 
slightly, while for holes a difference in pulse height and length is observed. 



The errors given in Tables [T] and [2] have been deter- 
mined by Minuit but should be seen as only approximate 
due to the high non-linearity of the system. 

To evaluate the influence of systematic uncertainties in 
the thickness and depletion voltage the same set of mea- 
surements was fitted, assuming ± 1 % difference in thick- 
ness or ± 5 % difference in depletion voltage. The dif- 
ference in the resulting mobility parameters is within the 
range of the errors calculated by Minuit. 

The penetration depth of 660 nm light changes by about 
50 % in the investigated temperature range. Compari- 
son of simulations with the longest and shortest penetra- 
tion depth and otherwise unchanged parameters showed 
no significant change in the overall pulse shape, thus the 
simulation neglected the temperature dependence of the 
penetration depth. 



Optical photons of 660 nm wavelength (1.87 eV) pro- 
duce so called hot charge carriers as the band gap of silicon 
is approximately 1.12 eV at room temperature. The simu- 
lation assumes charge carriers in thermal equilibrium with 
the crystal lattice. This is considered to be a good approx- 
imation as the thermalization of the hot carriers is usually 
very fast compared to the pulse duration. 

6.1. Fitted mobility parameters 

The fitted values of the ohmic mobilities at 300 K re- 
produce the accepted values [l7[ within the quoted errors. 
For the temperature coefficients reported values range from 
-2.3 to -2.6 for electrons and -2.2 to -2.8 for holes [17]. The 
best fit value of holes lies within the quoted interval, while 
the electron value of this work lies just outside the quoted 
interval (rounded to the second digit it would be -2.3 and 
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Figure 5: Comparison of different parameterizations of the saturation 
velocity proposed by other authors. In the pubhcation of Jacoboni et 
al. ^ two parameterizations are given. The one parameterizing the 
electron saturation velocity directly has been labeled '2nd model'. 



thus at the edge of the quoted interval), although it is 
lower than the value of -2.0 which has been reported by 

i- 

Saturation velocity parameterizations as function of 
temperature have been provided by many authors jsl, lisl, 
ll9|, however the values differ by about 20 %. The val- 
ues for electrons of this work are in the range of values 
proposed by others, while the values for holes are higher 
than reported by [3|, [isl, ll9|, however values as high as 
0.96x10^ cm/s at 300 K have been found by Smith and 
Frey [20], which is even larger than the value determined 
in this work. A comparison is shown in Figure [5l 

Except from the value of [3] for <111> orientation, no 
parameterizations of /3 values as function of temperature 
and crystal orientation are available. 

A comparison of mobilities and resulting drift veloci- 
ties determined in this work with those proposed by other 
authors is shown in Figure [6l A distinctive difference in 
the transition behavior between ohmic mobility and satu- 
ration velocity between the different parameterizations is 
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Figure 6: Comparison of different parameterizations of the mobility 
and drift velocities at T = 300 K. Electrons are marked by the full 
symbols, holes by the open symbols. Vertical black lines indicate the 
approximate minimum and maximum field encountered in the mea- 
surements of this work. Both electrons and holes show an anisotropy. 
Comparing <111> to <100> orientation the electron drift velocity 
and mobility decreases, while the hole drift velocity and mobility 
increases. 



observed. Especially for electrons a lower mobility than 
predicted in [3, 5] is observed (for both crystal directions). 
In the transition region hole mobilities in the <111> di- 
rection tend towards the values proposed by [5] , while the 
mobilities in the <100> direction tend towards the values 
proposed by [3|. In the saturation region the wide spread 
in the proposed velocities (shown in Figure [5j) dominates, 
especially at lower temperatures (not shown). 

7. Discussion and Conclusions 

A new parameterization of the field and temperature 
dependence of electron and hole mobilities in bulk silicon 
was presented for <111> and <100> crystal orientation. 

The differences between the values derived from the pa- 
rameterization presented here and the parameterization of 
Jacoboni et al. [3[ are approximately ± 15% in the inves- 
tigated temperature regime. The differences as function of 
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Figure 7: Difference between mobilities determined in this work and 
those of ^ at T = 300 K. Vertical black lines indicate the approxi- 
mate minimum and maximum field encountered in the measurements 
of this work. Maximum deviations at this temperature are zb 12%. 



electric field are shown in Figure [71 While the parameter- 
izations for electrons in the <111> show an almost con- 
stant difference of- 5%, the other parameterizations show 
distinctive differences. The differences to the parameteri- 
zations proposed by others [j, |5|, |6|, |7| are even larger and 
are not presented here. 

No samples of <110> orientation were available, thus 
no parameters for this direction could be determined. The 
electron parameters for this orientation provided by Cas- 
toldi and Rehak [21] (/i^ = 1394cmVVs, <^^ = 1.85x10^ cm/s 
^<iio> _ 2.51) differ greatly from the parameters pre- 
sented here. However they have been obtained with very 
low fields (up to 700 V/cm), in which case the simple pa- 
rameterization used in this work might need some correc- 
tions. 

The presented results resolve the ambiguities in the 
temperature dependent mobility parameterizations for the 
<111> orientation j3|, y, |5|, |7| and present a parameteri- 
zation for the <100> orientation. The results obtained 
enable more precise simulations of the charge collection in 
silicon sensors. The study is part of the investigation of 
the plasma effect, which occurs when the density of the 
electron hole pairs produced by the radiation exceeds the 
bulk doping density of the silicon. Its understanding is of 
relevance for the design of silicon pixel sensors for experi- 
ments at x-ray free-electron-lasers. It turned out, that the 
precision of the literature data on the transport parame- 
ters was insufficient to describe the reference data at low 
charge carrier densities. Therefore more precise data was 
required for a quantitative simulation of the plasma effect 
Hi, 23]. 
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